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Abstract 



As a toy model for dynamics in nonequilibrium quantum field the- 
ory we consider the abelian Higgs model in 1 + 1 dimensions with 
fermions. In the approximate dynamical equations, inhomogeneous 
classical (mean) Bose fields are coupled to quantized fermion fields, 
which are treated with a mode function expansion. The effective equa- 
tions of motion imply e.g. Coulomb scattering, due to the inhomoge- 
neous gauge field. The equations are solved numerically. We define 
time dependent fermion particle numbers with the help of the single- 
time Wigner function and study particle production starting from in- 
homogeneous initial conditions. The particle numbers are compared 
with the Fermi-Dirac distribution parametrized by a time dependent 
temperature and chemical potential. We find that the fermions ap- 
proximately thermalize locally in time. 
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1 Introduction 

One issue that is of considerable interest in nonequilibrium quantum field theory 
is equilibration and thermalization. An understanding of this is crucial for the 
knowledge of e.g. the rate at which a quark-gluon plasma forms in heavy ion 
collisions, or the rate at which a Bose gas thermalizes locally in time during 
evaporative cooling. In cosmology, the study of (pre)heating of the universe at 
the end of inflation has become a topic of its own. 

Real time evolution of quantum fields can in general not be solved exactly. 
A popular approach in nonequilibrium field theory is then to use large N or 
Hartree-like approximations, in which coupled equations for mean fields and their 
fluctuations are solved self-consistently [1], ||, 0, f|, |5|, || 0j. In many treatments 
available in the literature, the mean fields are taken to be homogeneous. Ther- 
malization may be difficult to achieve in these approximations [fj, || £|, [5|, |(| (see 
in this respect also ||). Another approach is to treat the dynamics of the low 
momentum modes classically || [10| . 

In this paper we study numerically the real time dynamics in the abelian 
Higgs model in 1 + 1 dimensions, extended with fermions. The choice of model 



is motivated by electroweak baryogenesis [!1|, e.g. according to the scenario in 



ref. [fL2|. We use a large N approximation in which the Bose fields are treated 



as mean fields and the fermion fields play the role of fluctuations. Previous 
numerical studies of fermions in real time have been restricted to homogeneous 
mean fields, with 0, |6| or without [[U| back reaction. Instead, in our case the 
mean fields are inhomogeneous and they obey the full non-linear classical field 
equations including the back reaction of the fermions. This back reaction is 
represented by mode functions which obey the Dirac equation in the presence of 
the mean fields. Loosely speaking, the effective equations describe a collection of 
quantum mechanical particles (the fermion modes) coupled to classical fields (the 
mean fields). The particles therefore interact and scatter amongst themselves as 
in classical electrodynamics. In particular Coulomb forces, screened by the Higgs 
field, are expected to play an important role. 

We focus on the following dynamical problem: initially, all the energy is 
contained in a few long wave lengths of the Bose fields, and the fermions are in a 
vacuum state. This is a nonequilibrium situation and in course of time energy is 
transferred towards the fermionic degrees of freedom, producing fermion particles. 
A motivation is partly given by inflation, where at the end of inflation the energy 
stored in the inflaton field is transferred to other degrees of freedom, leading to 
the (pre)heating of the universe (see e.g. |3], || [13], PH). 



In the presence of spacetime dependent mean fields the identification of par- 
ticles is usually based on the concept of adiabatic particle numbers [Q, |3|, [|, [)]. |(| . 
If the model can be described in terms of weakly coupled quasiparticles, it should 
also be possible to deduce a particle number directly from correlation functions. 
In a gauge theory, simple correlation functions such as the fermion two-point func- 



tion are not gauge invariant and the calculation would have to be performed in a 
fixed gauge. This introduces an ambiguity whether the resulting particle number 
depends on the gauge. However, physical observables calculated in terms of these 
particle numbers should come out gauge independent. Another possibility is to 
use a gauge invariant modification of the two-point function. There are several 
ways in which this can be done and a simple one is by including a parallel trans- 
porter between the two fermion fields at different points in space. This gives a 
close connection to the Wigner function. 

The Wigner function has a long history in transport theory and quantum 
kinetic theory [|I5 |. It has been used to find approximative ways to deal with 



dynamical issues in out-of-equilibrium plasmas. Recently, the Wigner function 
has been used quite extensively to derive transport equations that contain the 



same physics as the Hard Thermal Loops in 3 + 1 dimensional gauge plasmas [16]. 
A comparison between so-called covariant and single-time formulations is given 
in 



17|j . In this paper we use the single-time Wigner function, not for dealing 
with the dynamics but merely as a tool for identifying the particle number. We 
solve the microscopic dynamics numerically in the large N approximation and 
then calculate the single-time Wigner function. After averaging over space and a 
short interval in time to smoothen the effect of oscillating Bose fields, we extract 
the particle number by comparing the Wigner function with that of a free fermion 
gas. This approach is obviously only valid when the theory is relatively weakly 
coupled and a description in terms of quasiparticles makes sense. 

The paper is organized as follows. In section |2| we introduce the model and 
give the effective equations of motion. A discussion of the particle number from 
the fermion two-point function and the Wigner function is given in section ^. In 
section 01 we present numerical results for the real time evolution and the time 
dependent particle number, starting from a nonequilibrium initial state. To see 
whether the particles are effectively thermalized, we compare the particle numbers 
with the Fermi-Dirac distribution depending on time dependent temperature and 
chemical potential. The results are summarized in the Conclusion. The equations 
are solved numerically on a lattice in space and time. Some aspects of the lattice 
implementation are briefly described in appendix |A]. Details relating the Wigner 
function and the gauge fixed two-point function are given in appendix ||. 

We refer to our previous paper [0] for more details on the model, the effective 
equations of motion, and a discussion of the lattice implementation. 

2 Action and mean field equations of motion 

We consider the abelian Higgs model in 1 + 1 dimensions, with rif = 2 flavours 
of fermions. The action reads 

S = -Jd 2 x ^F, u F^+(D^yD^ + \(4 > *4 > -v 2 /2) 2 



+^7 M («9 M - iqA^i + - X)(G^fCV*^ - GMWi 



Here i 7 )^ = d fJ- A l/ — d v A^, D^ = <9 M — iA M , C denotes the charge conjugation 
matrix and ip = ip^/3. An explicit representation for the gamma-matrices is 7 1 = 
01, C — (3 — ry° = cr 2 , and 75 = — 7°7 1 = 03. Space is a circle with circumference 
L and the Bose (fermion) fields obey (anti)periodic boundary conditions. Other 
conventions can be found in . The choice of this model, in particular the form 
of the Yukawa coupling, is motivated by electroweak baryogenesis, but we will 
not elaborate on those aspects here (see instead 0). The fermion field carries a 
flavour index i = 1,2 = nj and the Yukawa coupling G, depends on the flavour: 

*i, ipi -> e iq k 
so that the Yukawa term is gauge invariant. In this paper 
we restrict ourselves, however, to massless fermions and put the Yukawa term to 
zero, G — 0. 

The global symmetry ipi — > e luJ15 ipi is broken in the quantum theory due to 
the axial anomaly, and 



G\ = —G 2 = G. Local gauge symmetry acts as <fi — > e^0, ipi — > e %qK ip, 



An + d^, with q 



1 

2' 



d»J5 



f A-K ^ v 



n f d,C^ 



with the axial current j£ = iipi'j^'j^ipi and the Chern-Simons current C M = 
€ ixu A u /2it. A change in the axial charge is proportional to a change in the Chern- 



Q 5 (t)-Q 5 (0) 

with Q 5 (t) = / dxj®(x,t), and C{t) 



n f (C(t)-C(0)), 
-- -(27V)- 1 J dxA 1 (x,t). 



Simons number 

(2.1) 

As in the elec- 
troweak theory, integer values of the Chern-Simons number label the non-trivial 
vacua of the bosonic theory. The q = | charged fermions change the periodicity 
of the classical bosonic ground state: only vacua with Chern-Simons numbers 
that differ an even integer (instead of any integer) are connected by large gauge 
transformations. The vacua are separated by finite energy barriers, the sphaleron 
configuration at half integer C. 

As an approximation to solve the dynamics, we treat the Bose fields as mean 
fields and the fermion fields as fluctuations. A formal derivation can be given by 
duplicating the fermion fields N times and taking the large N limit, after a proper 
rescaling of the coupling constants and fields [I], |7]]. The resulting equations are 
(in the A = gauge) 



9 2 A(x,t) =e 2 [j 1 h (x,t) + (j 1 f (x,t)) 
d%(f)(x,t) = Dl<f>(x,t) - 2\ [\<p(x,t)\ 2 - v 2 /2 

together with Gauss' law 

d 1 d A 1 {x,t) = -e 2 \j° h (x,t) + (j%x,t)) 



(f>(x,t), 



(2.2) 
(2.3) 

(2.4) 



The scalar contribution to the current, j£, is given by 

j£{x, t) = i(D»<p(x, t))*(p(x, t) - i(j>*(x, £)Z>'0(x, t). (2.5) 

The fields A\ and <fi represent mean fields, which can be inhomogeneous in general, 

A\{x, t) = (A 1 (x, t)), (j)(x, t) = ((j)(x, t)). 

The fermion contribution to the current is 

j?{x,t) = fy$ i (x,t),'fil> i {x,t)] i 

and its expectation value represents the fermion back reaction. The commutator 
ensures that the current is odd under charge conjugation. 

The most demanding part of this approach is the calculation of the back 
reaction. In the case of homogeneous mean fields, the usual manner is to use a 
mode function expansion. We use the same method in the inhomogeneous case 
as well. 

The back reaction is calculated by expanding the fermion fields in a complete 
set of eigenspinors of the initial Dirac hamiltonian at t — 0. This leads to a mode 
function expansion, 

i/>i(x, t) = Y1 [bpiU pi (x, t) + d ] pi v pi (x, t) . (2.6) 

p 

The coefficients b pi , d pi are time independent and represent the annihilation (cre- 
ation) operators of particles resp. antiparticles at the initial time. They obey the 
usual anticommutation relations 

{ bpi 7 bp'i'} = 5 PP > 5 hi , {d pi , d p >i> } = 5 pp i 5u> , 

and zero for others. The expectation values of the creation and annihilation 
operators in (|2.6j ) determines the initial state for the fermions. We take vacuum 
initial conditions for the fermions, so that the only non-zero expectation values 
are 

Mji) = (dpidl) = 1. (2.7) 

The fermion contribution to the current can now be expressed in terms of the 
mode functions, 

{jf(x, £)) = - J2 9 fapifo fyl^Upiix, t) - Vpi(x, t)i^v pi {x, t)\ . (2.8) 

The time dependence is carried by the spinor mode functions u and v, which are 
solutions of the Dirac equation in the presence of the gauge field, 

id u pi (x,t) =H D [A 1 (x,t)]u P i(x,t), (2.9) 



and similar for v P i. The Dirac hamiltonian has the usual form 

H D = -ia 1 ^ - iqA x ), a 1 = -7V. (2.10) 

The initial conditions for the mode functions are given by 

UpiyX, U ) y— 6 Upi, VpiyX, V) /~T ^pii 

V L yL 

with Upi(v p i) positive (negative) energy spinors of the Dirac hamiltonian at t = 0: 
p labels the momentum eigenstates of the initial Dirac hamiltonian. It takes the 
values 

27T 1 

due to the antiperiodic boundary conditions on the fermion fields. For these 
initial conditions with vacuum expectation values fl2.7|), the current ( |2.8| ) vanishes 
at t = 0. 

If the mean fields are restricted to be homogeneous, plane waves e ±ipx mul- 
tiplying the mode functions can be factored out at all times, and each mode 
function is the solution of an ordinary differential equation. In the more general 
case that the mean fields are allowed to be inhomogeneous, this is not possible, 
and a large set of partial differential equations has to be solved. In that case the 
label p looses its interpretation for t > 0, and it simply labels a complete set. 

The equations of motion are free from ultraviolet divergences and the sum 
over the mode functions in ( j2.8p is finite. In the higher dimensional case charge 
renormalization would be necessary [[!]]. For a study of renormalization in the 
presence of the Yukawa term we refer to our previous paper 0. The energy of 
the system is the sum of the energy in the Bose fields and the fermion fields. 
The energy of the fermions at time t is given by the expectation value Ep(t) = 
^2 x {ipi(x : t)7iDil)i(x, t)) and can be expressed as a sum over mode functions. This 
sum is quadratically divergent. In the numerical implementation this is regulated 
by the lattice cutoff: the effective equations are solved using a discretization on 
a lattice in space and real time (see appendix 0). We renormalize the energy by 
subtracting the bare fermion energy at t — 0. 

3 Particle number and the Wigner function 

Observables which may give insight about what happens on the microscopic level 
during time evolution are current, charge and energy densities. These are gauge 
invariant quantities whose physical meaning is clear and unambiguous and they 
are straightforward to calculate. In addition to this, it has been shown in the 
case of homogeneous mean fields that more insight can be found with the concept 
of adiabatic or instantaneous particle number. 



The adiabatic particle number is based on an expansion of the (interacting) 
quantum field at each instant in time in terms of a complete set of eigenf unctions 
of the Dirac hamiltonian at time t (we suppress flavour indices here) 

ij}{x, t) = E [bE(t)u E (x, t) + d E {t)v E {x, t)] . (3.1) 

E 

The coefficients of these eigenf unctions are then identified with time dependent 
creation and annihilation operators. In ( [3.1|) the eigenf unctions are denoted with 
u E ,v E , and the complete set is labeled with E. A Bogoliubov transformation 
relates the expansions (|2.6| ) and (|3.1| ). The instantaneous particle numbers at 



time t are then defined by the expectation values 

N E (t) = (b E (t)b E (t)), N E (t) = (d E {t)d E {t)). 

In the case of homogeneous mean fields the instantaneous eigenf unctions can be 
calculated analytically, in terms of plane waves. For the inhomogeneous case the 
diagonalization will in general have to be performed numerically, which makes 
the calculation rather involved. Furthermore, the instantaneous eigenmodes of 
the Dirac hamiltonian are not easily interpreted in terms of familiar concepts 
like quasiparticles, e.g. concerning the relation between E and the quasiparti- 
cle momentum. Therefore we now discuss another approach, using correlation 
functions.^ 

The identification of particles can be motivated in general (for a weakly cou- 
pled system) by comparison with the noninteracting case, and for free fields the 
two-point function provides all the available information. Fields in equilibrium 
are homogeneous in space and time. In the interacting situation, one can look for 
subsets of the system which are effectively in equilibrium. Such subsets will only 
exist locally in space and time, and one way to enforce homogeneity is averaging 
over small subset-volumes of spacetime. This then leads to a definition of particle 
numbers locally in space and time. The optimal size of the subsets is not known 
a priori, and this may introduce a certain ambiguity. The dependence on the size 
of the subsets can be studied by changing it. 

One way to implement this idea is by studying the equal-time two-point func- 
tion 

S{x,y,t) = &{x,t)${y,t)), (3.2) 

in particular in the following manner 
S av (p,t) = — f t+tm 2 dt' S{p,t'), S(p,t) = fdze- ipz j fdxS{x,x + z;t). 

t-av Jt—tetv/2 J Li ■> 

(3.3) 



1 Correlations functions have been used in other contexts as well: see e.g. |L8J for a study 
of defect formation during nonequilibrium evolution with correlation functions. 



The time average over a short time interval t av will smoothen oscillations in the 
Bose fields, and instead of averaging over smaller spatial subsets, we choose here 
to integrate over the complete spatial volume.0 Hence, the particle numbers 
will be defined locally in time only. This simplifies the analysis and has the 
positive effect of diminishing the fluctuations. As a side remark, it can be checked 
explicitly that this approach reduces to a commonly used definition of particle 
number for the case of homogeneous mean fields |5|. || . 

The fermion correlator is not gauge invariant, but the concept of quasipar- 
ticles is not gauge invariant in the first place, although using them to compute 
observables such as energy and pressure should give gauge invariant answers. 
Therefore, before averaging over space and time we need to fix the gauge. A 
suitable choice is the temporal Coulomb gauge, A = 0, diA% = 0, in which 
the remaining freedom of 'large' gauge transformations is removed by requiring 
C = -LA 1 /(2tt)6(-1,1]. 

There are of course other possibilities. One is to use gauge invariant modifi- 
cations of the fermion correlator (|3.2|), for example, 



W(x,y;t) = (i>(x,t)i>(y,t)U(y,x;t)), (3.4) 

where U is the parallel transporter 



U(y, x; t) = exp 



-iq / dsAi(s,t) 
Jy 



The path is along the shortest straight line connecting x and y. Other paths are 
possible as well, which leads to a class of gauge invariant modifications of the 
two-point function ( |3.2| ). The arbitrariness in choosing a gauge is now replaced 
by choosing a gauge invariant correlator. The above W(x, y\ t) is closely related 
to the single-time Wigner function |T5| , I7 |. In terms of the relative coordinate 



z = x — y and the center-of-mass coordinate X = l(x+y), the single-time Wigner 
function is given by a spatial Fourier transform of Q3.4p with respect to z, 



W(X,p;t) = Jdze~ tpz W(X + ^z,X- -z;t). 
We shall use its spatial and temporal average, 

W m {p,t) = ±- [ t+tav/2 dt'W(p,t'), W(p,t) = \ fdXW(X,p;t), (3.5) 

£av Jt—tav/2 h J 

similar as for the two-point function, ( |3.3|) . In appendix [B] we show that the 
correlators S(p,t) and W(p, t) are actually closely related by mere interpolation 

2 In equilibrium the two-point function in ( |3.3| ) does not depend on x at all because of 
translational invariance. Slightly out of equilibrium, the long wave length inhomogeneities 
around equilibrium in the two-point function (or rather in the Wigner function, see below) can 
be used to derive transport or kinetic equations by applying a gradient expansion |15[ |16[ . 

8 



between the discrete momenta p. In the following we shall continue with the 
Wigner function. 

In the spinor decomposition (in 1 + 1 dimensions), four functions appear 

W(X, p; t) = T{X, p; t) + i^V, (X, p; t) + z 75 P(X, p; t) , (3.6) 

that are real, due to the elementary property W\X,p;t) = /3W (X , p; t) (3 . We 
shall use the obvious notation T aw {p, t), etc., for temporally and spatially averaged 
coefficients. 

For fermions characterized by momentum p and mass m, the free field expres- 
sion for the Wigner function is given by 

W irec (X,p) = {1-N P - N- P ) m ~^ lp + iz 7 °(l - N p + AL p ), (3.7) 



E p = Jm 2 +p 2 . 

The arbitrary (anti)particle occupation numbers are given by N p (N p ). Since 
this is independent of X, we now take as a possible definition of time dependent 
particle numbers N p (t), N p (t) and effective mass m p (t), the solution of the three 
equations 

1 1 r i m (t) 

J- av (p, t) = -Tr W av (p, t) = - 1 - N p (t) - N„ p (tj\ -, »\ } , (3.8) 



2 ""^'^ 2 



^ml(t)+p 2 



V lav (p,t) = JWVMP,*) = ^ fl " ^p(*) - ^-p(*)1 / = . .(3-9) 



2 -'""«-'-' 2 



^Jml(t)+p 2 



VoavGM) = ^Trz 7 °^ av (p,t) = i [l - iV p (t) + N. p (t)] . (3.10) 

The explicit solution of these equations is given in appendix JTJ. We have added 
a subscript p on the effective mass parameter to indicate a possible momentum 
dependence. The interpretation of m p as a mass is of course best when m p (t) is 
momentum independent, as in the free field case. A comparison with free Wigner 
function ( |3.7| ) gives zero for the fourth function V av (p,t) = — Tr27 5 W av (p, t). 
Hence, this coefficient cannot be used to find new information on the particle 
numbers or effective mass, but it serves as a consistency check. 

Note that when the mean fields are homogeneous, particles and antiparticles 
with momentum p can only be produced in pairs, and N p = N_ p . In general, 
they can be different. It should be clear that we use the Wigner function only 
to identify the instantaneous particle number at time t (i.e. after averaging over 
the short interval t av ). Therefore we consider the single-time Wigner function 
instead of the covariant one. We do not attempt to solve the dynamics using the 
Wigner function approach. 

Our goal is now to solve the effective equations for the Bose fields and the 
spinor mode functions. The Wigner function can be expressed in terms of the 



mode functions using the expansion (|2.6| ), and from a calculation of JF av (p,i) 



Vna V (p,t) and V av (p,t) the particle number can be extracted. This will be the 
subject of the following section. 

4 Particle production 



We solve the closed set of effective equations ( |2.2h ]2~5l, [2~8| - |2.10| ) numerically. In 



this paper we take the following parameters: the dimensionful parameters are 
related by A/e 2 = 0.25, eL = 3.2 and the dimensionless parameter in the scalar 
potential is v 2 = 8. Then the volume is not very large but also not too small: in 
terms of the tree level boson masses, m^L = y2\vL « 6.4, tuaL = evL « 9.1. 
Furthermore, the couplings are fairly weak: e 2 /m 2 = 0.25. Some details on the 
lattice implementation and choice of parameters can be found in appendix |A[ 

The inhomogeneous Bose fields are initialized such that all the energy initially 
resides in the long wave lengths. We take vacuum field configurations <p(x, 0) = 
v/y2, Ai(x, 0) = 0, and use space dependent momenta. The results presented 
below are obtained starting with 

do4>(x, 0) = 3ecos — — \- 2zecos — — , d Ai(x, 0) = Gauss' law + e 2 . 

The electric field d§A\{x, 0) is determined almost completely by the Gauss law 
( |2.4|) , up to a constant. 

The Chern-Simons number C and the axial charge per flavour (Qs) /n/, which 
agree in accordance with the anomaly equation (|2.1| ), are shown during time 
evolution in fig. |l[ Time is given in units of 1/e, which has the dimension of 
mass in 1 + 1 dimensions. We see one sphaleron transition. The energy of the 
Bose fields, the fermion fields and the total conserved energy are shown in fig. 
0. Initially, the (renormalized) fermion energy is zero, because of the vacuum 
initial conditions for the fermions. We see that there is energy transfer towards 
the fermion degrees of freedom. 

To analyse how the energy is distributed in the fermion subsystem we have 
calculated the four functions !F(p, t), V M (p, t) and V(p,t) which appear in the 
spinor decomposition of the Wigner function (|3.5| , |3.6| ),p| and averaged them over 
a relatively short time interval et av = 4. This interval is comparable with the 
elementary periods present in the Bose system: et^ = 2ixe/rriA w 2.2, et^ ~ 3.1. 
Below, the time et refers to the center of the interval et av . With regard to 
the spatial momentum dependence, p is written in terms of the dimensionless 
combination ap, where a is the lattice spacing. Note that the continuum regime 
is where the momentum p is small with respect to the cutoff ir/a, e.g. \ap\ £ 0.5 
(see appendix |A|). 

3 The results for the particle numbers presented below are normalized for one flavour, i.e. 
divided by rif = 2. 

10 



We solve the time dependent particle number from the three equations that 
are suggested by the free Wigner function (|3.8H3.ltiD , in the combination 



N 

1 v avp 



(t) = ~ (N P (t) + JV^(t)) . 



The result for the averaged particle number is shown in fig. || The particle 
numbers of the modes with \ap\ > 1 are consistent with zero, up to numerical 
precision. These modes are not excited at all. As can be seen from fig. [|, essen- 
tially only the physical (i.e. without lattice artefacts) fermions (with \ap\ ^ 0.5) 
are excited. Furthermore the distributions are much smoother than those ob- 
tained with homogeneous mean fields || |(J. Note that the distribution functions 
are not symmetric under p — > —p. This will be discussed below. 

With respect to the time dependent mass m p {t) in (|3.8H3.10|) , we find no 
significant deviation from the free dispersion relation E = \p\Q Finally, we find 
that Voav(p, t) — 1/2 and V av (p, t) are consistent with zero for momenta \ap\ > 0.5. 
For smaller momenta they fluctuate around zero. The amplitude in the case of 
Va,v(p,t) is smaller than 0.005. 

In order to see whether the fermion subsystem effectively thermalizes locally 
in time we compare it with the Fermi-Dirac distribution function, depending on 
a time dependent temperature T(t) = l//3(t) and chemical potential fJ,(t), 

fp{t) = exp[P(t)(E p -^t)q 5p )] + l (41) 

The chemical potential is coupled to the axial charge density, and we use the 
notation q 5p = sign(p) = sign(7 5 ) to denote the chirality Note that this term 
breaks the symmetry p — > —p. In thermal and chemical equilibrium N p (t) = 
N- P (t) = iV avp (£) = f p (t). 

To compare the particle numbers with the Fermi-Dirac distribution, we plot 
log(A r a ^ 1 — 1) versus ap. In equilibrium this would result in a straight line 
/3(t)[|p| ± n(t)]. The results are shown in fig. |[ We see that the points seem 
to lie approximately on straight lines, with small deviations present. The lines 
through the points are obtained from a straight line fit. Note that the data 
denoted with et = 2 represent the averaged particle number in the first inter- 
val, < et < 4, which contains less than 1.5 oscillations of the Chern-Simons 
number (see fig. [I]). Already here the low momentum modes show approximate 
thermalization. 

From least square fits to the data we can find the time dependent tempera- 
ture and chemical potential, with errors. The effective temperature is obtained 
separately from p > and p < distribution functions and is shown in fig. [5[ 
The effective chemical potential is shown in fig. |[ In this case we present the 
average of the chemical potential obtained from the p > and p < distribution 



4 In the analysis we use of course the lattice dispersion relation. 

11 



functions, since the fluctuations are rather large. The errors correspond to the 
least squares fit. 

If the fermions are in thermal and chemical equilibrium, with possibly slowly 
varying time dependent temperature and chemical potential, and interactions 
are neglected, the relation between (Q5), Ep, T and /i is as follows. For massless 
fermions, E p = \p\, we find that in that case the axial charge is given by 



(Q 5 ) r°° dp 



nfl ^ 



1 



e 0(E p -fiq5 p ) _|_ I e /3(E p +nq 5p ) _|_ J 



71 



With the help of the anomaly equation, (Q 5 ) = rifC, this gives /x in terms of 
C. Likewise, the renormalized fermion energy at finite temperature and chemical 
potential can be calculated, 

E F fTC 2 1 2 \ TCC 



L »4«r + S"J- " = -L- (42) 

Given C and Ep, we can predict T and \x assuming the fermions are in equilibrium 
and neglecting interactions. This predictions are shown in figs. |5| and |6] as well. 
(Here the error is the standard deviation of the average in an interval et av = 4, 
treating the data as uncorrected. ) We see that the equilibrium temperature 
obtained from ( |4.2j ) lies significantly below the effective temperatures obtained 
from the particle distribution functions. Presumably this reflects the fact that 
the fermions are of course not really free. 

The chemical potential from the fits has large fluctuations and errors. This is 
possibly due to finite size effects. The time independent equilibrium result (|4.2| ) 
suggests that the chemical potential will be sensitive to time dependence of the 
Chern-Simons number (see fig. [I]). This dependence should decrease when the 
volume L gets larger. 

5 Conclusion 

We presented numerical results for fermion particle production in the presence of 
inhomogeneous time dependent Bose fields. The quantum field dynamics was ap- 
proximated by effective mean field equations of motion for the Bose fields coupled 
to quantized fermion fields, represented by mode functions. The particle num- 
ber was extracted from the single-time Wigner function, using the free Wigner 
function as a guideline. For strongly coupled theories, such an approach would 
presumably not be possible. 

A comparison of the particle number with a time dependent Fermi-Dirac dis- 
tribution, i.e. containing a time dependent temperature and chemical potential, 
showed that the produced particles are approximately in thermal equilibrium, 
already after the first few oscillations of the Bose fields. The effective temper- 
ature increases when more energy is transferred towards the fermion degrees of 
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freedom. Note that these results are quite different from those obtained with only 
homogeneous Bose fields, in which case particle distributions are typically non- 
thermal. Due to the inhomogeneous gauge field, the equations studied here con- 
tain classical (screened) Coulomb scattering of the particle-like modes, whereas 
this is absent in the case of homogeneous mean fields, for which mode coupling 
is severely reduced. 

In this paper our main focus was on fermion particle production. An analysis 
of the complete system, e.g. of the time dependent energy distribution of the Bose 
fields, remains to be done. 

This study may be relevant for inflationary scenarios. Preheating of the uni- 
verse is usually analysed by coupling the homogeneous inflaton field to the parti- 
cle mode functions, leading to resonance band structures and highly non-thermal 
particle distribution functions (see e.g. references given in the Introduction). The 
stability of these resonance bands and the time scales of thermalization will be 
affected when inhomogeneities are taken into account. 
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A Lattice implementation 



We solve the effective equations numerically, using a formulation of the theory on 
a lattice in space and time. For details concerning especially the fermion fields 
we refer to 0, here we repeat only the necessary. We use a lattice in space with 
N sites and lattice spacing a, such that L = Na. The index p labeling the mode 
functions then takes a finite number of values 

p=^E(n-i), ne{-±N + l,...,±N}, (A.1) 

and the total number of mode functions is given by 2nfN = 4N. On the lattice 
the dispersion relation is modified and reads for massless Wilson fermions J?| 



E p = JSp + m,p, s p = a 1 sinap, m p = a x (l — cosap), (A. 2) 

where m p is the Wilson mass'. The maximal fermion momentum on the lattice 
equals ap = tt(1 — 1/N) ~ n, in units of the lattice spacing. For such high 
momenta lattice artefacts are important. The physical modes are those with 
momenta that are small in lattice units, since then E p = \p\(l + 0(a 2 p 2 )). A 
rough guideline is \ap\ < 0.5, for ap = 0.5 the correction term is approximately 
1%. The time step is denoted with oo- The lattice parameters are N = 64 and 
a / a = 0.005. 

13 



B Numerical calculation of the Wigner function 

In this appendix we describe how we calculate the Wigner function numerically 
on the lattice. We also compare the single-time Wigner function with the gauge 
fixed equal-time two-point function. 

While the definition of the Wigner function is manifestly gauge invariant, it 
is particularly convenient to calculate it in the (completely gauge fixed) Coulomb 
gauge. We define this (suppressing the time dependence) by 

A\ = A l {x) + di9 c (x) 7 d x A\ = 0, -1< C c < 1, (B.l) 

where the gauge fixed Chern-Simons number is C c = —LA\/(2 , k). This com- 
pletely fixes the gauge. In particular the freedom under large gauge transforma- 
tions is fixed by the last requirement in ( |B.1| ). We recall that due the q — ~ 
charged fermions, C — 1 is not gauge equivalent to C = 0, but e.g. to C — — 1. 
The fermion fields transform accordingly ip{x) — > ip c {x) = exp[iq8 c (x)]tp(x). The 
parallel transporter in the Wigner function is in this gauge simply U(y, x;t) — 
exp [iqA\(y — x)]. 

The Wigner function, averaged over the lattice in space, becomes 

W(p, t) = \ £ £ e-^ + ^>(MX + \z, t)MX - ±z, t)). (B.2) 

Up to the factor exp [—iqA\z\, this is identical to what would have been obtained 
when started from the Coulomb gauge fixed propagator. To evaluate it, we first 
calculate the Fourier transform of the gauge fixed two-point function 

S(p, t) = £ e-*"i Y.(MX + \z, t)UX - \z, t)), (B.3) 

which is the discrete Fourier transform of an antiperiodic function. Hence we 
know S(p,t) at a discrete set of p values, namely those given by ( |A.1| ). Using 
inter- and extrapolation, the value of S(p, t) at other values of p can be found 
as well. It follows from comparing ( |B.2| ) and ( |B.3| ) that the Wigner function is 
given by W(p, t) = S(p + qA\, t). 

The gauge fixed propagator shows discontinuous behaviour when the Chern- 
Simons number crosses the boundaries set by fixing the freedom under large 
gauge transformations in ( |B.1| ). We expect these discontinuities to become less 
important when the volume is increased. On the other hand, the presence of the 
time dependent phase factor in ( |B.2| ) leads to a smooth behaviour of the Wigner 
function when the Chern-Simons number changes during time evolution. 

Finally, the explicit expressions for the time dependent particle number N p (t), 
the antiparticle number N p (t) and the effective mass m p (t) in terms of the func- 
tions J-gcipjt) and Vna, v (p,t) can be obtained as follows. First note from 
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and (OR that 



1 /. „ .-, \2 



KM t) + VUp, t) = ^(l-N p - N. PJ , 
which gives |1 — N p — N- p \. The sign can be found from (|3.9|) , and 



iV avp (t) = ^(N p + N_ p ) = sign[pV lav (p, t)] yjFl(p, t) + Vl v (p, t) + i 

The individual (anti)particle numbers (instead of the sum) can be found by 
adding or subtracting 

-Vo, v (P,t) + \ = \(N p -N_ p ). 

The effective mass, or more generally the dispersion relation, can then be found 
from 

provided Vi av (p, t) 7^ of course. However, Vi av (p, t) vanishes only when N &vp {t) — 
I equals zero as well, since the momentum p on the lattice is always different from 
zero, due to the antiperiodic boundary conditions in space (see (|A.1|) ). Note that 
since Tsv{p,t) and Vi av (p, t) are both real, J r a v(p,t)/Vi av (p,t) = —m p (t)/p is real 
as well, which implies that m p (t) is real. 

References 

[1] F. Cooper, S. Habib, Y. Kluger, E. Mottola, J. P. Paz and P. R. Anderson, 
Phys. Rev. D50 (1994) 2848. 

[2] Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper and E. Mottola, Phys. 
Rev. Lett. 67 (1991) 2427; F. Cooper, S. Habib, Y. Kluger and E. Mottola, 
Phys. Rev. D55 (1997) 6471. 

[3] D. Boyanovsky, H. J. de Vega, R. Holman and J. F. J. Salgado, Phys. Rev. 
D54 (1996) 7570; D. Boyanovsky, D. Cormier, H. de Vega, R. Holman, 
A. Singh and M. Srednicki, Phys. Rev. D56 (1997) 1939. 

[4] D. Boyanovsky, C. Destri, H. J. de Vega, R. Holman and J. Salgado, Phys. 
Rev. D57 (1998) 7388. 

[5] Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper and E. Mottola, Phys. 
Rev. D45 (1992) 4659. 



15 



[6] D. Boyanovsky, M. D'Attanasio, H. J. de Vega, R. Holman and D. S. Lee, 
Phys. Rev. D52 (1995) 6805; J. Baacke, K. Heitmann and C. Patzold, Phys. 
Rev. D58 (1998) 125013. 



[7; 

[8 

[9 

[10 

[11 
[12 

[13 

[14 
[15 
[16 

[17; 

[1 



G. Aarts and J. Smit, Nucl. Phys. B555 (1999) 355. 

L. M. A. Bettencourt and C. Wetterich, Phys. Lett. B430 (1998) 140. 

D. Y. Grigoriev and V. A. Rubakov, Nucl. Phys. B299 (1988) 67. 

S. Y. Khlebnikov and I. I. Tkachev, Phys. Rev. Lett. 77 (1996) 219. 

For recent reviews, see e.g., V. A. Rubakov and M. E. Shaposhnikov, Usp. 
Fiz. Nauk 166 (1996) 493; M. Trodden, |hep-ph/9803479| . 



J. Garcia-Bellido, D. Grigoriev, A. Kusenko and M. Shaposhnikov, 
Eiep-ph/9902449 . 



P. B. Greene and L. Kofman, Phys. Lett. B448 (1999) 6; G. F. Giudice, 
M. Peloso, A. Riotto and I. Tkachev, JEEP 08 (1999) 014. 

L. Kofman, A. Linde and A. A. Starobinsky, Phys. Rev. D56 (1997) 3258. 

D. Vasak, M. Gyulassy and H. T. Elze, Ann. Phys. 173 (1987) 462. 

J. P. Blaizot and E. Iancu, Nucl. Phys. B390 (1993) 589; Phys. Rev. Lett. 
70 (1993) 3376; Nucl. Phys. B417 (1994) 608. 

P. Zhuang and U. Heinz, Ann. Phys. 245 (1996) 311; S. Ochs and U. Heinz, 
Ann. Phys. 266 (1998) 351. 



D. Ibaceta and E. Calzetta, [hep-ph/98 10301 . 



16 



_>-, 







20.0 30.0 

et 



50.0 



Figure 1: Chern-Simons number C and axial charge per flavour (Qs)/nf versus 
et. The lines fall on top of each other, as expected from the anomaly equation. 
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Figure 2: Energy in units of e in the Bose fields (Eb), fermion fields (Ep) and 
their conserved sum Et versus et. 
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Figure 3: Time dependent particle number iV av versus ap, the momentum in 
lattice units, at three different times. The particle numbers are obtained from a 
time average over an interval of length et av = 4 and et denotes the center of the 
interval. 
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Figure 4: Transformed particle number log(Ai r av 1 — 1) for the data of fig. || In 
equilibrium this would give /3(£)[|p| ±/x(t)]. The lines are obtained from a straight 
line fit. 
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Figure 5: Time dependent effective temperature versus et, obtained from least 
square fits (with errors) from the transformed particle numbers as in fig. |], for 
p > 0,p < separately. The line denoted with Ep, C is obtained from ( |4.2| ), i.e. 
assuming that the fermions are in thermal and chemical equilibrium. 
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Figure 6: Time dependent effective chemical potential versus et, obtained as 
in the previous figure. Here the result is averaged over the chemical potentials 
obtained from the p > and p < fits. The line denoted with C is again obtained 
from ([Oft. 
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